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ABSTRACT 



A model is developed which relates the structure of a 
combat unit to its expected effectiveness during a time inter- 
val. Sensitivity to variations of structure is incorporated 
into the model by adapting effectiveness concepts from reli- 
ability theory, and by constructing the attrition process as 
a continuous parameter Markov chain. The model is highly 
sensitive to the judgment of the user, due to parametric 
flexibility both in the attrition process, and effectiveness 
measurement. Two solution methods are given. Examples are 
employed to illustrate the model's use, and proposals for 
refinement are offered. 



2 



TABLE OF CONTENTS 



I . INTRODUCTION 5 

II. RELIABILITY CONCEPTS 9 

III. MEASUREMENT OF EFFECTIVENESS 19 

IV. UNIT PERFORMANCE MODEL 38 

V. CONCLUSION 57 

BIBLIOGRAPHY 60 

INITIAL DISTRIBUTION LIST 61 

FORM DD 1473 62 



3 



LIST OF FIGURES 



1. Series System 14 

2. Parallel System 15 

3. Two-Out-of-Three System 16 

4. Parallel System with Irrelevant Component 17 

5. Example Outcome of Actual Effectiveness Random 

Variable and Value Random Variable 27 

6. Coherent Systems Associated v;ith the Absolute 

Effectiveness Function of Example 1 35 

7. Constant Matrix A 52 

8. Functional Matrix [si - B] 53 



4 



I. 



INTRODUCTION 



The purpose of any operational model is to aid a decision- 
maker by relating his possible alternatives to their consequent 
results, thus increasing the knowledge upon which his decision 
is based. Therefore a model, in addition to validly represent- 
ing the process or system involved, must be capable of dis- 
playing the contrast, in terms of meaningful criteria, of even 
subtle variations of choice. Sensitivity to variables outside 
the area of the current decision is not particularly valuable, 
and in fact may produce detrimental complication. In short, a 
model should ideally be a reflection of the problem of its user. 
It is for this reason that distinct decision problems often 
require the use of different models of a single system. 

This dependence of model upon problem is nowhere more 
apparent than in the study of combat dynamics. As an example 
a single battle situation may be modeled using game theory if 
the decision involves strategy, or using search theory if 
alternative fire distributions are being studied, or using 
Lanchester theory if weapon characteristics or tactical move- 
ments are being considered. Each m.odeling approach has been 
found useful for certain problem types but useless for others. 

The purpose of this thesis is to propose an alternate 
modeling viewpoint for those problems in which the structure 
of a combat unit must be related to its effectiveness. The 
term structure is meant to denote the personnel and equipment 
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of the unit, the roles of these elements, as well as the 
channels employed to receive and transmit information and 
instructions. As an example, a decision-maker who must 
choose for field test a small number of prototype tank killer 
teams from all the infinite organizational possibilities, 
faces a structural decision problem, and might profit from 
the model to be proposed. 

The importance of structure derives from the fact that 
v;hen a leader or communicator becomes a casualty, his removal 
from the structure produces tv/in transformations in both the 
effectiveness and vulnerability of the unit. These transforma- 
tions are determined at least in part by the casualty's posi- 
tion in the organization. As an illustration consider a rifle 
squad in defense whose squad leader has just become the first 
casualty. His function is assumed by a fire team leader, v;ho 
must nov; receive and transmit orders, direct the fire of his 
squad and search his sector for targets. The effectiveness of 
the squad is reduced not only because one fev;er weapon is 
firing but because a less highly trained leader is in charge. 
The vulnerability has increased for the team leader because he 
now must increase his exposure in order to direct the whole 
squad and observe its whole sector. The vulnerability of the 
whole squad is perhaps increased because it is being directed 
with less expertise. 

Thus, to know the condition of a unit it is not enough to 
know the number of casualties without knowing v;ho they are. 
Perhaps the strongest testimony to this viewpoint is the 
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universal military principle that enemy leaders and radio 
operators are the targets of highest priority in any type of 
engagement, and the simultaneous effort to conceal visual 
signs of their identity. If the validity of this approach 
has been justified the next step is to outline the model 
requirements . 

The decision-maker v;ill be assumed to have available to 
him a facility for military judgment which will include his 
own assessment of unit and individual vulnerability, and his 
own criteria for measuring unit effectiveness on the basis of 
its surviving structure. The mission of the model will be to 

1) provide a framework in which the decision-maker's 
assessment of threat can be specifically quantified, 

2) translate the threat assessment into a stochastic 
description of the behavior of the candidate 
organization, 

3) provide a framework for expressing the decision- 
maker's effectiveness criteria as an algebraic 
function of the unit condition, 

4) derive a stochastic description of effectiveness of 
the candidate organization, using the effectiveness 
function . 

Fortunately, the complete process of model construction is 
unnecessary, since there already exists an established field 
of study which treats the concept of structure. The field is 
stochastic reliability theory, which has historically devel- 
oped in answer to the need for relating the design of complex 
systems of mechanical and electronic components to system 
effectiveness. The comparison of a system of components to 
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a unit of fighting men may seem rash. However, through a 
process of adaptation, a model v/ill be derived which may be 
considered more consistent while still taking advantage of 
reliability concepts. 

Additional reliance on the theory of continuous parameter 
Markov processes will permit treatment of significant inter- 
dependencies in the attrition process. 

The resultant model is one with capacity for immediate 
application, as well as promising opportunities for improve- 
ment and refinement. 
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II . RELIABILITY CONCEPTS 

Reliability theory deals v;ith collections of inter- 
dependent components. These collections are called systems. 

A component is any entity v/hich can be described at any time 
as either functioning or failed, depending upon its ability 
to perform its assigned mission. Components are generally 
assumed to be in one of these two conditions. No partially 
failed conditions are considered. 

The system v/hich is discussed in this thesis is an engaged 
combat unit whose components are included in one or more of 
the following categories: 

1) personnel, 

2) weapons and equipm.ent, 

3) channels of information and corrim.unication . 

The number of components in the unit is represented by 
the variable n called the order of the unit. Each component 
is identified with an integer in the set {l,...,n}, called its 
index . 

The performance random variable of the ith component is 
defined as follows: 

X. (t) = 0 if the ith component is in the 
failed condition at tim.e t, 

=1 if the ith component is in the 

functioning condition at time t. 

The component reliability of the ith component is defined 
as follows: 
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r^(t) = P[X^(t) = 1] . 

Since there are n components of interest, it is an obvious 
extension to define the performance random vector as follov/s: 

X(t) = (Xj^(t) , . . . ,Xj^(t) ) . 

The outcome space of X(t) is the set of all binary n-vectors. 

Call this set . This set has N members where N = 2^. 

n 

It is often necessary to partition the outcome space 
into m subsets, a,,..., a , where 1 < m < N, and define a 
partition function by assigning a unique real number b^ to 
all vectors in each subset a^ , so that 

a (x) = b^ if and only if x e a^ , 

Each vector is a member of exactly one of the subsets. 

The first partition function of interest is the state 
function s (x) , which is constructed by partitioning into 
N subsets, each containing one vector. The value assigned to 
each one-member subset is one of the integers 1,...,N. 
Therefore 

s (x) = j if and only if x e a^ and b^ = j . 

It is assumed here that in all cases 

s (1) =1 and 

s (0^) = N where 

_1 = 1) and 

0 = ( 0 , . . . , 0 ) . 
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The composition of the state function with the perfor- 
mance random vector results in definition of the state random 
variable 

S ( t) = s (X (t) ) . 

The common practice in reliability theory is to consider 
the entire system to be either in the functioning or failed 
condition at any time. Predictably the system's condition is 
described by a binary random variable. However, v;hen the sys- 
tem is a combat unit, such a measure is too coarse to be of 
value. Still, since this binary concept is of value in the 
adaptation process, it is introduced belov/. 

The structure function 52$ (x) of a system of order n is a 
partition function constructed by forming two subsets: 

1) a, , the set of all vectors x in Q , such that if 
X(t) = X, then the system is functioning, 

2) a2 , the set of all vectors y. such that if 

X(t) = y, the system, is in the failed condition. 

Consistent with the use of binary indicators, 

^(x) = l(=bj^) if and only if x e a^^ 

= 0(=b2) if and only if x e a2* 

Composing the structure function with the performance 
random vector results in definition of the structure random 
variable 

<^(t) = 0 {X(t) ) . 
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The system reliability is defined to be 



R(t) = P[$(t) = 1] . 

If more than one system is being considered then it is neces- 
sary to place a subscript on the structure function, struc- 
ture random variable and system reliability. 

A system is called coherent if its structure function 
satisfies three conditions: 

1) (z5(0) = 0, 

2 ) ( 2 ^( 1 ) = 1 , 

3) If X £ then 0(x) <. . 

Essentially a coherent system, is one which fails if all its 
components fail, functions if all its components function, and 
whose condition is not impaired by improving the condition of 
any of its components. 

A component is called irrelevant if its condition never 
has any effect on the value of the structure function. 

Components are assumed to be of two varieties. Renev;able 
components alternate between the functioning and failed condi- 
tions. For example, tactical radios are communication channels 
whose functioning is interrupted by atmospheric disturbance or 
jamming. Nonrenewable components never return to the function- 
ing condition after entering the failed condition. Personnel 
are nonrenewable components of a unit. 

A common graphic method of illustrating the relationship 
of component performance to system performance is the use of 
the block diagram. Examples of these diagrams are shown in 
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Figures 1-4. The method is derived from the common practice 
of representing electrical circuits. Each block corresponds 
to the like indexed component of the system. A potential 
flow is imagined, attempting to pass through the diagram 
along the connecting arcs. The flow is able to pass through 
a component block only if the component is functioning. If 
the flov; is able to traverse the diagram over any path then 
the system is functioning and its structure function equals 
one. Irrelevant components are represented by broken blocks 
totally unconnected to the system. 

Figure 1 shows the block diagram for the system of three 
components which functions if and only if all three components 
are functioning. This is called a series system. Figure 2 
shows a system which fails if and only if all its components 
fail. This is called a parallel system. Figure 3 shows a sys- 
tem which functions if and only if at least tv/o of the compo- 
nents are functioning. This is called a two-out-of-three 
system. As in this case, a single component may be represented 
by more than one block, with the understanding that all blocks 
with the same index are always in the same condition. Figure 
4 shows a system with components one and two in parallel, with 
component three irrelevant. It should be rem.erabered that the 
block diagram shows only structural relationship, and implies 
nothing concerning actual location. 

A final concept, not unique to reliability, but useful 
here, is that of Laplace transformation. Given a function 
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Figure 1. Series System 
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Figure 2. 
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Figure 3. Two-Out-of-Three System 
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Parallel System v;ith Irrelevant Component 
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c(t), the Laplace transform of c(t). is defined as follo^vs: 



L {c (t) } = c (s) 



O 

/ 



exp (-St) 



c(t) dt, 



a function of s. The function c(t) is called the inverse of 
the Laplace transform c(s). 
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III. MEASUREKENT OF EFFECTIVENESS 



In order to make any description of unit performance 
useful to a decision-maker, it is necessary to translate that 
description into some measure of effectiveness or goal attain- 
ment. Such a measure is useful only if it is consistent v/ith 
the decision-maker's standards of value, and is in a form, 
sufficiently simple to permit reasonable comparison of compet- 
ing alternatives. The method developed here, of constructing 
such a m.easure begins with the designation by the user of a 
relatively simple instantaneous indicator function based upon 
his intuitive value judgments. This function then serves as 
a basis for derivation of a one dimensional deterministic 
function, which measures unit quality over a variable time 
interval . 

The actual effectiveness of an engaged unit at a point in 
time is a very complex function of numerous variables, includ- 
ing friendly and enemy force levels, and many immeasurable 
physical and psychological factors. The simplifying assumption 
is made here that effectiveness depends only upon the vector of 
component performance and the length of time elapsed since the 
beginning of the engagement. The enem.y performance is not 
considered directly, but the designation of an effectiveness 
function should reflect consideration of a hostile enemy envi- 
ronment, based on a previously determined scenario. The 
friendly unit is always assumed to have all components func- 
tioning at time zero. 
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The decision-maker's first step is to define a partition 



function on The set is partitioned into subsets a^ , 

i=l,...,m/ where vectors in the same subset are judged to 
represent equal effectiveness unit conditions. Vectors in 
different subsets represent outcomes which are not equally 
effective. This partition is performed without reference to 
time. To each subset a^ is assigned a dimensionless real 
number g^ , called its absolute effectiveness level. The 
function thus defined is called the absolute effectiveness 
function g(x) where 

g(x) = g^ if and only if x e a^, 
and the partition is constructed so that 

It is assumed without serious limitation that 



(3-1) 




g(0) =0, 


(3-2) 






(3-3) 


if X < ^ 


then g (x) < g (y_) 



These assumptions merely require that a completely destroyed 
unit be considered completely ineffective, that a completely 
undamaged unit be considered maximally effective, and that 
loss of a component never improve effectiveness. 

Since the vector menibers of and the integers in the 
set {1 ,...,n} are placed in one-to-one correspondence by the 
state function s (x) , it is possible to designate a companion 
function to g(x) called the alternate effectiveness function 
defined on {1,...,N} such that 
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if and only if 



d(j) = 

It follows that 



s (x) = j and £ ^ ^ • 



g (x) = d (s (x) ) . 

This function is useful in later calculations. 

The construction of the actual effectiveness function is 
of extreme importance to the successful employment of the 
model. The user should take special care to assure that the 
absolute effectiveness levels truly represent his judgment, 
not only v;ith respect to ordinal relationships, but to ratio 
relationships as well. 

A unit's real effectiveness depends upon hov; long it has 
been engaged. One factor to be considered is the change in 
enemy situation due to losses or reinforcement as provided in 
the scenario. Another consideration is that in missions of 
fixed length, such as delaying actions, casualties are less 
damaging, the later they occur. To allow for such trends, the 
user's next step is to define a time variance function , f(t), 
and a resultant actual effectiveness function, 

h (x, t) = g (x) f ( t) . 

If in a specific application the variance effect is not con- 
sidered, then 

f(t) =1, t >L 0 

and therefore 

h (x, t) = g (x) , t ^ 0 . 
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The function of g (x) can be expressed in another equiva- 
lent manner which shov/s the close relationship with reli- 
ability theory. Esary, Proschan, and Walkup [Ref. 3] point 
out that a function such as g(x), having a finite range can 
be written as 

(3-4) g(x) = (g^-gj^) (x) + 

+ (g -g 1 ) 0 1 (x) 

^m ^m-1 ’^m-1 — 

where the functions (x) are defined such that 

(3-5) (is) “ if only if g (x) <. g^ 

(3-6) (is) ~ i if only if g(x) > g^. 

It is now shown that each of the functions jz5^(x) is the struc- 
ture function of a unique binary coherent system composed of 
the components of the unit, in which at least one of the 
components is not irrelevant.' Furthermore the systems are 
related in such a way that no system is functioning if any 
lower indexed system has failed, and no system is failed if 
any higher indexed system is functioning. 

Theorem ^ 

The functions in (3-4) are structure functions of 

coherent systems of the n components of the parent system. 
Proof 

It must be shov/n that for i=l,...,m, 

1 ) 9i^{0) = 0 , 

2 ) ^^( 1 ) = 1 , 

3) f f — Z. then 0^ (x) <. (^) . 
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The proof is by contradiction. The above three points are 
each in turn assumed false, and a result, contradicting one 



of the assumptions made in defining g (x) , is shown to follow. 

1) Assume that (0_) 0 for some i. Then it must be 

that (0_) = 1, but then from (3-4) 
g(0) 1 ^ ° 

contradicting (3-1). Thus (0_) = 0/ for i = 1, . . . ,m-l 

2) Assume that 0^^ (^) 1 for some i. Then it must be 

that 0^(^) = 0, but then from (3-4) 

g(l) < (gj-gi) + ... + = g^,. 

g (1) < g 

^ — ^m 

contradicting (3-2). Thus 0j^(^) = 1, for i=l,...,m-l. 

3) Assume that there exist vectors x and y in such 

that X <. ^ and that for some i, 0j_(x) > • Since 

the function is binary, 0^_(x) = 1, and (x.) =0* By 
(3-6), 0^(x) = 1 implies that g(x) > g^. By 3-5), 

~ ^ implies that 9 — '?i* Therefore 

g (y) < g(x) and ^ — Y. 

contradicting (3-3) . Thus if x £ y , then 
0j_(x) £0^(y), for i = 1,..., m-l.|| 

Theorem ^ 

The functions 0^(x), i = 1,..., m-1, are related as follows 
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1) if 0}^(x) = 1/ then ei j (x) = 1, for j = 1,..., k-1, 

2) if 0j^(x) 0/ then fi-M = 0, for j = k + 1,..., m-1. 



Proof 

The 

1 ) 



proof is again by contradiction. 

Let ^j^(x) = 1 some vector x. If k=l then the 

result is trivially true. If k > 1, assume that 
(z5.(x) = 0 for some j < k. By (3-6), 0, (x) = 1 implies 

3 “ K — 

that g(x) > g, . By (3-5) , . (x) = 0 implies that 

— K 3 — 

g(x) <. g • . Since g. < g, , it follov/s that 
— D 3 

g(x) < gj^ < g(x) , 

an obvious contradiction. Therefore 0.(x) = 1 for 

3 - 



j = 1 , . . . ,k-l . 



2) Let 0j^(x) ■= 0 for some vector x. If k = m-1, 
result holds trivially. If k < (m-1), assume 
0.(x) = 1 for some j > k. By (3-5), jzJ, (x) = 

J — K — 

that g(x) <. g, . By (3-6) , 0. (x) = 1 implies 

— K 3 — 

g(x) > g.. Since g. > g^ , it follov/s that 
g(x) < g^ < g(x) , 

an obvious contradiction. Therefore (x) = 

j = k + 1, . . . , m-1 . II 
Theorem. 3_ 

For any vector x, g(x) = g^^/ if and only if ^j^(x) 

^i-l(x) = 1. 



the 

that 

0 implies 
that 



0 , for 



= 0 , and 
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Proof 



1) Let g(x) = for some x. Then by (3-5), g(x) = 

that 0^(x) = 0. By (3-6), g(x) = 
g^^ > implies that = 1. 

2) Let 0^(x) = 0, and (x) = 1 for some x. By (3-5) 

(x) = 0 implies that g(x) <. g^. By (3-6) , 

0j^_2^(x) = 1 implies that g(x) > Therefore it 

follows that 

9i_i < g(>i) ^ 

but the only value in the range of the function v/hich 
satisfies the inequality is g^. Therefore g(x) = g^.H 
The ultimate importance of this method of expressing g(x) 
is that the study of a very complicated multivalued absolute 
effectiveness function can be simplified to the study of a 
linear combination of simple dependent binary structure func- 
tions of the type commonly treated in classic reliability. 

It follows that the actual effectiveness function can be 
written as follows: 

h(x,t) = f(t) { (g 2 -gj^) {z5j^(x) + ... + 

Thus far the functions introduced are deterministic with 
domain Since is also the outcome space of the perfor- 

mance random vector, X(t) , each of the functions can be 
composed with this vector to define analogous random variables. 
The absolute effectiveness random variable is defined by the 
relation 

G(t) = g(X(t) ) , 
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( 



or equivalently 




where 



^>^(t) = 0^(X(t)), i = 1,..., m-1. 

The actual effectiveness random variable is defined by the 
relation 



These effectiveness random variables have the advantage 
that they are defined through a process which assures their 
consistency v;ith the decision-maker's standards. Their pri- 
mary limitation is that they represent only instantaneous 
effectiveness. That is for every nonnegative value of t 
there is defined a distinct random variable H(t). In Figure 
5, the bold horizontal line depicts a possible curve of out- 
comes of H(t), where f(t) is assumed to equal one for all non- 
negative t. 

Since combiat engagements are extended over indefinite time 
periods, it is desirable to aggregate the effectiveness curve 
into a variable which measures accumulated effectiveness over 
an interval. To fill this need the value random variable is 
defined by the relation 



H(t) = f (t)G(t) . 



t 



t 



V(t) 




H(s) 




G (s) f (s) ds . 



o 



o 



In Figure 5, the value random variable is represented by the 
shaded area below the curve of H(t). 
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Figure 5. Example Outcome of Actual Effectiveness Random 
Variable and Value Random Variable 
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The value random variable is a measure of the type 
desired, but in general it will have a complicated distribu- 
tion, too unwieldy to be of practical use in decision making. 
A relatively simple measure of effectiveness, which still 
preserves many advantages of the value random variable, is 
the quality function defined by the relation 

q(t) = E[V(t) ] , 

the mean or expected value of the value random variable. 

To compute the quality function in a specific case, note 

that 



q(t) = E 



H(u) du 



E[H(u)] du 



o . 
t 



f (u) E [G (u) ] du. 



where the interchange of expectation and integration in the 
second line is justified by application of Fubini's Theorem, 
since H (u) is a family of nonnegative random variables. The 
validity of this interchange is established by Royden [Ref. 



12 ] . 

There are two approaches for calculating q(t). The first 
arises from the fact that 
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E[G(u)] = ^ g(x) P[X(u) = X |x(0) = 1] 



X 



N 



where . 



^d(i) p^^(u), 

i=l 



(u) = P[S(u) = i I S(0) = 1] 



is called the interval probability function from state one to 
state i. Thus the first result is that 



q(t) = 



t 

r 




- N 


“ 




1 f (u) 

) 




^d(i) 




du 


o 




L i=l 






N 


p 


t 




-i 


E^(i) 




r f (u) 


p^^(u) du 


i=l 


_ 


y 

O 







N 



(3-7) = ^d(i) I(i,t) 



i=l 



where 



(3-8) I(i,t) = f(u) P 3 ^j^(u) du. 
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The second approach follows from relation (3-4), 

m-1 

E[G(u)] = E[ ^ + ^.(u)] 

3 = 1 ^ 



m-1 

= ^ E[.f.(un 

3 = 1 

m-1 

3=1 

m-1 

“ Z <93+1-33 1 '’ 3 <''l' 

3 = 1 

where (u) is the system reliability of the jth associated 
system, computed on the basis of all components functioning 
at time zero. And the second result is that 



q(t) = 





- m-1 


J f (u) 


Z < 93 + 1 '« 3 ’ 


0 


- j=l -1 


m-1 


r t 



du 



“ Z <^ 3 +r^:' / 



f (u) Rj (u) du 



j = l 
m-1 



o 



(3-9) = ^ (g^^^-g^) J(j,t) 

j = l 



where 



(3-10) 



J(j,t) = / f(u) Rj (u) du. 



Now since the functions f(t) and g(x) or d(i) are assumed 
already defined by the decision-maker, all that is required for 
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calculation of the quality function is either set of 
functions {p^^(u),..., {R^(u),..., The 

determination of the functions is treated in Section IV. 

The quality function is the primary measure of effec- 
tiveness proposed here. It is a one dimensional determ.inistic 
function, yet it is based upon consideration of the random- 
ness of unit performance. Since it is derived from functions 
defined by the decision-maker, it is consistent with his 
implicit judgment. And since it is a function of time, it 
displays sensitivity to the length of engagement. 

In order to illustrate the use of the quality function in 
predicting unit effectiveness, this section v;ill be closed v;ith 
an example. 

Example 1 

Consider a tank killer team consisting of the following 
components : 

COMPONENT INDEX 

Team Leader 1 

Gunner 2 

Recoilless Rifle 3 

It is assumed that the unit is engaged in a large scale battle 
as part of a large force. The performance random vector for 
this unit, X(t) , has an outcome space consisting of the 

eight binary 3-vectors. Let the state function be defined as 
follows ; 
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s{ (1,1,1) } = 1 
s{ (1,0,1) } = 2 
s{ (0,1,1) } = 3 
s{ (1,1,0) } = 4 



s{ (1,0,0) } = 5 
s{ (0,1,0) } = 6 
s{ (0,0,1) } = 7 
s{ (0,0,0) } = 8 



The decision-maker is now faced with the problem of 
constructing the absolute effectiveness function g(x). Sup- 
pose the hypothetical team is considered to have two missions. 
At any time the v;eapon is functioning, the primary m.ission is 
to detect enemy tanks and bring them under fire. If the 
weapon is in the failed condition, the secondary mission is to 
detect enemy tanks and direct the fire of adjacent teams 
against them. From experience, field test, or some other 
information source, the user determines that when both person- 
nel and the weapon are intact the team is capable of bringing 
effective fire on eight tanks, per hour. He decides to base 
his absolute effectiveness function on this fire capability, 
and his first step is to set g{ (1,1,1)} = 8. He feels that 
in most units the leader and gunner will have equal ability in 
firing the weapon. Therefore v.'hen either member becomes a 
casualty v;hen the weapon is functioning, the decrease in effec- 
tiveness is about the same. The user estimates the reduced 
capability to be about five tanks per hour, and defines 
g{ (0,1,1)} = g{ (1,0,1)} = 5. When the weapon fails, the team 
performs the secondary fire direction mission. The user deter- 
mines that both members together can detect and bring fire on 
as many as three tanks per hour. However, if the leader alone 
survives, he has a capability of tv;o tanks per hour, while the 
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gunner alone has a capability of only one per hour because he 
has less experience and training than the leader in directing 
fire. In the two cases v/here both members are casualties the 
capability is, of course, zero. Thus the decision-maker 
defines his absolute effectiveness as follows: 



g{ (1,1,1) } = d(l) = = 8 

g{ (1,0,1) } = d(2) = = 5 

g{(0,l,l) } = d(3) = g^ = 5 

g{ (1,1,0) } = d(4) = g^ = 3 

g{ (1,0,0) } = d(5) = g 3 = 2 

g{ (0,1,0) } = d(6) = g 2 = 1 

g{(0,0,l)} = d(7) = g^ = 0 

g{ (0,0,0) } = d(8) = g^ = 0 



The only requirement to be met in practice is that the 
resulting function faithfully reflect the true values of the 
decision-maker. Additionally, when several units are being 
compared, a separate function must be defined for each one, 
and however they are constructed, it must be the case that 
states of different units adjudged by the user to be equally 
effective have the same value of g(x). 

As indicated in (3-4) , g(x) can be expressed as follows: 
g(x) = (1-0) ^^(x) + (2-1) 02 (x) + (3-2) ^ 3 (x) + (5-3) (x) + 

(8-5) 0^(x) where the structure functions 0^ (x) are defined 
in the table below. 
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X 5^1 ^2^-^ 

( 1 , 1 , 1 ) 1 1 

( 1 , 0 , 1 ) 1 1 

( 0 , 1 , 1 ) 1 1 

( 1 , 1 , 0 ) 1 1 

( 1 , 0 , 0 ) 1 1 

( 0 , 1 , 0 ) 1 0 

( 0 , 0 , 1 ) 0 0 

( 0 , 0 , 0 ) 0 0 



^3(x) ^4(x) 

Ill 
110 
110 
10 0 
0 0 0 

0 0 0 

0 0 0 

0 0 0 



The block diagrams of coherent systems associated with these 
structure functions are shown in Figure 6. 

The decisicn-m.aker also wishes to define a time variance 
function for the tank killer team which enhances the actual 
effectiveness of every state the longer the engagement lasts 
He feels that enemy losses and confusion make friendly fire- 
power more valuable later in the battle. He consequently 
defines 

f(t) = 2-e“^. 



Thus f(0) = 1 and at the outset h(x, 0) = 9 (^i^ / after one 

hour h(x, 1) = (1.633)g(x). 

Assume the following conditional distribution for S(t) is 
given . 



P 

P 

P 

P 



ll(t) 


-It 

= e 


-3t -4t 

Pis(t) = e -e 


12 


-6t -It 
= e -e 


/4-^ _ -2t ^-3t 

Pj^g(t) e e 


13 


-5t -6t 

= e -e 


(t) = e"^-e 


14 


-4t -5t 

= e -e 


Pig(t) = l-e"'*^ 
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Coherent Systems Associated with the Absolute 
Effectiveness Function of Example 1. 
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By (3-8), 



1(1, t) 
1(2, t) 
1(3, t) 
1(4, t) 
1(5, t) 
1(6, t) 
1(7, t) 
1(8, t) 



9 



168 

3 
70 

1 

15 

7 

60 

1 

4 

5 

"6 
2t - 



7 

3 

2 



5 

1 



2 

2 



-7t 




1 


-8t 










8 


e 






-6t 




3 


^-7t 


1 


-8t 




+ 


7 


-e - - 


8 


e 


-5t 




1 


"6t 


1 


-7t 




+ 




- e 


8 


• e 


-4t 




3 




1 


-6t 




+ 


5 


- e “ 


~T~ 


• e 


-3t 


+ 


3 


-4t 

- e 


1 


-5t 
• e 



, ,,-3t 

e + e 



1 -4t 

— e 



^ 3 -2t 

2e + — e 



+ 3e 



2 

-t 



1 -3t 

-^5— e 



1 -2t 

-?r— e 



And by (3-7) , 
8 



q(t) = d(i) I(i,t) 



i=l 

= 8 1(1, t) + 5 1(2, t) + 5 1(3, t) 
+ 3 1(4, t) + 2 1(5, t) + 1 1(6, t) 
+0 I (7, t) + 0 I (8,t) 



653 , 3 -8t 

280 8 ® 



6 -7t ^ 1 -6t 

T- e + -3- e 



5 ® 



1 -4t 

-r— e 



^ ^-2t 

— e - e 



Now the given interval probabilities are equivalent to 
the follov/ing reliability functions of the associated systems 
in Figure 6: 
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Rj^(t) = e 
R2(t) = 

R2(t) = 

R^(t) = 

-2t 

R^(t) = e 

By (3-10) , 



J(l,t) = — 


2 

7 




+ 


1 

8 


00 

1 

CD 


J(2,t) = 3 ’ - 


2 

5 


(D 

1 

ft 


+ 


1 

6 




J(3,t) = ^ 


1 

2 


-4t 

‘ e 


-f 


1 

5 


-5t 
• e 


J(4,t) = i|- - 


2 

3 


-3t 
■ e 




1 

4 




J(5,t) = 




- 2 t 

e 




1 

3 


-3t 
‘ e 



And by ( 3-9 ) , 

5 

q(t) = J(j,t) 

j = l 



= 3 J(l,t) + 2 J(2,t) + 1 J(3,t) 

+ 1 J(4,t) + 1 J(5,t) 



653 ^ 3 - 8 t 6 - 7 t 1 -5t 

280 7® "^3® 






/ 



which agrees with the 



previously computed result. || 
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IV. UNIT PERFORMANCE MODEL 



The quality function, q(t), has already been introduced 
as a measure of effectiveness. The next step is the introduc- 
tion of a model to translate the user's judgment and intuition 
about combat dynamics into a description of the stochastic 
performance process of the proposed unit structure. The pri- 
mary requirements for the model are that it be sensitive to 
structural alternatives, and that it provide results in a 
form suitable for input to the quality function. There are 
two alternative forms of performance information which can be 
used in calculation of q(t). If expression (3-9) is used, the 
inputs required are the reliability functions of the associ- 
ated coherent systems. If the calculation is performed accord- 
ing to expression (3-7) , the information needed is the set of 
interval probability functions {p^^ (t) , . . . (t) } . This 

second type of information is the type of performance descrip- 
tion which is most easily determined employing the model 
proposed in this section. 

The unit to be modeled is assumed to be a complex struc- 
ture which includes, as components, personnel, equipment, and 
channels of information. The assumption will be made that 
each of the components is either functioning or failed. At 
the opening of an engagement all components are considered to 
be functioning. Some components such as personnel are non- 
renev;able and do not return to the functioning condition once 
they leave it. Other components are renewable and can return 
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to the functioning condition repeatedly. For instance the 
information channel between a squad leader and a team leader 
which consists of voice communication can be interrupted by 
battle noise, then restored when the noise subsides. There- 
fore an engagement can be looked upon as a series of events 
in time which change the state of the unit, through attrition, 
interruption, restoration or some combination of these effects. 
The state of the unit remains unaltered between these events. 

One approach commonly used to portray this process is to 
consider the performance of each component to be probabilis- 
tically independent. This assumption is not acceptable for 
this application since it severely limits the consideration 
of dependence in the attrition process. 

In order to integrate structural sensitivity into the 
model, it is necessary that the interevent times and the na- 
ture of each event be considered to obey random distributions 
which depend upon the surviving structure. The following 
proposal for accomplishing this is derived in part from the 
work of Marshall and Olkin [Ref. 9]. 

The engaged unit is subject to a hostile environment v;hich 
consists of the delivery over time of one or more types of 
lethal projectiles, including nuclear v;eapons, artillery, and 
small arms fire. Two assumptions are made about these fires. 

1) The delivery of each type of projectile constitutes 

an independent time homogeneous Poisson process, 

having interarrival times distributed exponentially 

with expected value where k is the index of the 

k 

weapon type . 
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2) For each unit state and each weapon type the damage 
caused by each round obeys a probability mass func- 
tion: Pj^(j|i) = PtProjectile of type k causes 

transition to state j | unit is 
in state i] . 

From these assumptions it follows that when the unit is in 
state i, the distribution of arrivals of projectiles of type 
k, which send the unit into state j, is Poisson having expo- 
nential interarrival times with expected value — r-ri- . 

Now in general the unit in state i can be subject to 
transition into state j as a result of several threats. If 
there are z types of projectiles V7hich can cause the transi- 
tion, then the distribution of all such events is the super- 
position of z independent Poisson distributions, which is 
itself Poisson, having exponential interarrival times ^^7ith 

Pi CjT'il Ti' p;-(31T7t; • 

In addition to attrition events there are also interruption 
and restoration events which describe alterations in condition 
of the renewable components. It is assumed that for each state 
i all functioning renewable components have an exponential time 

to failure with mean where j is the state which the com- 

ponent ' s failure causes the unit to enter. Likewise, if state 

i has a failed renewable component, then is the mean of the 

°ik 

exponential time to restoration and consequent transition into 
state k. 
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Let be defined as the time interval from the unit's 

entry into state i until an event occurs causing transition 
into state j. is then the time until an attrition event, 

or an interruption-restoration event of the proper type 
occurs. This is again a superposition of independent Poisson 
processes resulting in a Poisson distribution having exponen- 
tial interarrival times with expected value. 



Pj^^C j li) T. 



+ 






Try 



+ 



0 . . 

ID 



1 

ID 



and survival function 

-(A . .) t 

P[T^j > t] = e . 

In general the unit in state i faces threats v;hich can 
cause transition to several different states. The outcome is 
determined by the type of the earliest state-changing event to 
occur. Let T^ by defined as ’the length of the time interval 
between transition a-1 and transition a, so that 



T^ min{T^^,T^2' • • • / 

if state i v;as entered at transition a-1. If transition from 
state i to state j is not possible, then T^^ is considered to 
be infinite, v/ith = 0. Howard [Ref. 6] refers to this 

type of model as a "competitive process" since it is equivalent 
to a process in which, upon entry into state i, the random 
variables T^^ are sampled and, in effect, compete for the 
earliest outcome. 
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The survival function of is derived as follows: 
PIT^ > t) = PtT.^ > t, T-2 > t T. > t) 



= > t! PIT.j > tl...P(T^^ > t) 



= e 



- dil + + • • • + ^ 



Thus T^ is distributed exponentially with mean. 



E[T^,] = 



a' 



r 



+ A . 






TT 



1 

r~ • 
1 



‘il i2 "iN 

Thus the original assumptions of Poisson arrival of projectiles, 
and exponential interruption-restoration times, imply that the 
times between state-changing events are also exponentially 



distributed with mean j— . 

i 

The validity of these assumptions must be considered in 
any application of the model. The assumption of Poisson inci- 
dence of fires is characteristic of relative stability in 
delivery rates. This distribution v;ould not be consistent 
with scenarios which included large-scale enemy reinforcement 
or temporary massing of fires. 

The exponential distribution of the interevent times, T^, 
implies the following result: 



-U^)(t+s) 

P[T > t + sIt > t] = — T = e ^ 

a ‘a ,• t 



= P [T > s] . 
a 

This result is referred to as the "memoryless" property or, in 
reliability theory, as the "no wear" property. The effect of 



42 



the property in this model is that a unit's probability of 
successfully remaining in its current state for one addi- 
tional minute, is completely independent of the amount of 
time already spent in the state. This property fails to ac- 
count for such factors as fatigue or depletion of supplies. 
However in any case where the effect of these factors is 
insignificant compared with the actual attrition and interrup- 
tion, the "no wear" property may prove to be acceptable. 

The major consequence of the exponential assumptions is 
that the Markov property holds. That is, for t ^ u ^ 0, 

P[S(t) = i|s(v) , u > V > 0] = P[S(t) = ilS(u)], 

or the distribution of future S(t) depends only upon the cur- 
rent state, and not at all upon earlier information. The 
disadvantage of this property is that it does not permit con- 
sideration of such phenomena as momentum and collapse of 
situation. In applications where such effects are considered 
important, alteration would be required. 

The fortunate consequence of the Markov property is that 
it qualifies the process, {S(t), t ^ 0}, as a continuous para- 
meter Markov chain, a type of process v;hich possesses a v;ell 
developed method of solution. 

The decision-maker using this model is presumed to have in 
mind a specific unit structure and a specific engagement 
scenario. From the scenario he determines the variety of 
projectile types to which the unit is vulnerable along with 
rates of fire Tj^. Then for each state i he must specify the 
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mass functions Pj^(j|i) for each type threat k. These mass 
functions are the primary ingredient of the model's sensi- 
tivity to structure. They should be based on such variables 
as each individual's level of training, the capacity in which 
he is serving due to attrition of key personnel, the degree 
to which he must expose himself to hostile fire as a result 
of his current position, and the current condition of channels 
of communication. The specification of the mass functions 
will be the most difficult part of parameterizing the model. 
But it is indispensible if worthwhile investigation of unit 
structure is desired. 

The fact that the Markov property holds leads to the 
derivation of two sets of simultaneous first order differen- 
tial equations v;hich can be solved for all the functions 
p^j (t) . These two sets of equations are known as Kolmogorov's 
forward equations and backward equations. The actual deriva- 
tion is not presented here but is given in Ross [Ref. 11]. 
Before presenting these equations, several notational concepts 
are introduced. 

There are N distinct functions Pj^^^ (t) , which appear in 
the functional m.atrix 



P(t) = 



Pll(t) ... 



p^l (t) . . . (t)_i , 
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where the function (t) appears in the ith row and jth 
column. For any t >. 0 the rov.’s of p(t) add to one. Since 



p. .(0) = 1 i=j 

= 0 ij<^ j , 

it follows that 



P(0) = I 



10 ... 0 

0 1 



1 0 

_0 . . . 0 ij . 



The functional matrix P'(t) is defined as follows 



P' (t) = 



P'll(t) ••• P'iM(t) 



IN 



P Nl ’** ^ 



NN 



where 



p'i.(t) = 



d p . . ( t) 
^1 j 

dt 



Using the parameters already defined in this section, the 
constant matrix A is defined as follows: 



A = 



21 



12 

0 



IN 



N1 



where is the reciprocal of the mean of and.X^j appears 

ith row and jth column. Since no transitions are considered 
from any state into itself X^^ = 0, for i = 1,..., N. 
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The constant matrix B is defined as follows: 



-Xi ^12 • • • ^IN 

^21 ” ^2 

• • • 
• • • 

^ N1 "^N 



N 

where X . = ) X . . , the reciprocal mean of the interevent time 

j=l when the unit is in state i. 

Kolmogorov's equations are stated in matrix form below. 

The backward equations are given by 



(4-1) 



P' (t) = B P(t) , 



with initial condition P(0) = I. The forward equations are 
given by 

(4-2) P' (t) = P(t) B, 

v;ith initial condition P(0) = I. The operation indicated on 
the right of each expression is matrix multiplication. The 
fact that B P(t) = P(t)B is of special interest since matrix 
multiplication is not, in general, comm.utative . The solution 
of the process will consist of solving the backward equations 
using Laplace transforms. 

The Laplace transform of any functional matrix A(t), is 
indicated by 

L{A(t) } = A(s) , 

which is simply the matrix of Laplace transforms displayed in 
the same order as in A(t). For any functional matrix A(t), 



46 








f 




(4-3) 



L{A'(t)} = A' (s) = sA(s) - A{0), 



an equation proved by Kreyszig [Ref. 8]. Also for any func- 
tional matrix A(t) and any constant matrix C, 

(4-4) L{C A(t) } = C A(s) 

(4-5) L{A(t)C) = A(s)C 

If Laplace transforms are taken of both sides of (4-1) 
the result is 

L{P' (t) } = L{B P (t) } , 
which by (4-3) and (4-4) is equivalent to 

sP (t) - P(0) = B. P (t) , 

or since P(0) = I, 

sP(t) - B P(t) = I, 
or 

[si - B] P(t) = I. 

Premultiplying both sides by the inverse of [si - B] gives, 

P(t) = [si - B]"^, 

which gives the solution for L{P(t)} in terms of the constant 
matrix B. The solution is complete when the transform is 
inverted to give P(t). 

Howard [Ref. 6] shows that 

L{P (t) } = [si - B]“^, 

if and only if 

P(t) = e®^ , 



47 




{ 



where 



I + Bt + |j + |y + 



Bt 

e 



I 



is called the matrix exponential function. This expression 
may be used to calculate approximations to the elements of 
P(t), but will not provide exact solutions. Such a solution 
can be arrived at by finding the inverse of [si - B] and 
using the method of partial fractions to translate the result 
into invertible Laplace transforms. 

Example £ 

Consider again the tank killer team of Example 1, with the 
same three components, all considered nonrenewable, and the 
same state function s (x) . 

The scenario of interest involves tV70 types of external 
threat. The first threat consists of small arms projectiles, 
delivered with a Poisson rate of fire = 50/hour. The 
second threat is artillery fire v;ith rate T2 = 10/hour. 

The decision-maker now determines the approximate proba- 
bility mass function for each state-threat combination. 

For state one and the small arms threat, the user first 
notes a single round cannot cause more than one casualty. He 
therefore defines P^(5 
since these transitions would entail multiple attrition. He 
feels that the gunner is about 50% more vulnerable to small 
arms than the leader because of his firing position, and the 
fact that the backblast of his weapon pinpoints his location. 
The weapon is considered only half as vulnerable as the leader 
because its small presented profile. The user thus estimates 
the following mass function: 
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P^(l|l) = .94 ° 

P^(2|l) = .03 Pj^Cejl) = 0 

Pj_(3|l) = .02 Pj^CVlD = 0 

P^(4|l) = .01 P^(8|l) = 0 



where the actual magnitudes may be based on field data, 
simulation, or intuition. 

For state one and the artillery threat multiple attrition 
is possible. The user feels that each artillery round has a 
capability to destroy all three components and defines 
P2(8|1) = .1 . The gunner is considered more vulnerable to 
artillery fire for identical reasons as cited under small arms. 
The decision-maker defines this mass function as follows: 



The probability of transition into state seven is considered 
to be zero, since a round placed well enough to eliminate both 
personnel would be expected to destroy the weapon as well. 
Transition to state six is similarly considered impossible. 

The actual magnitudes are again based on the user's judgment. 

Since there are no renewable components the parameters X ^ ^ 
can now be computed : 



p2(l|l) = .60 
P2(2|l) = .05 
P2(3|1) = .1 
P 2 (4 |1) = .05 



P2(5|1) = .1 
P2(6|1) = 0 

P2(7|1) = 0 

P2(8|1) = .1 
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(trivial transitions are ignored) 



^12 " Pi( 2 | 1 )t^ + 

= (.03) (50) + 
= 2 

^^2 ^ (3 1 1) + 

= (.02) (50) + 
= 2 

^14 “ Pj^(4|1)Tj^ + 

= ( .01) (50) + 
= 1 

(5 1 1) + 
= (0) (50) + (. 
= 1 

Xi6 = P^(6ll)x^ + 
= 0 

^17 = Pi(7ll)T3^ + 
= 0 

^18 "" Pj^(8|l)Xj^ + 

= 0 + (. 1 ) ( 10 ) 

= 1 



P 2 ( 2 | 1 )x 2 
(. 05 ) ( 10 ) 

P 2 ( 3 | 1 )x 2 
(. 1 ) ( 10 ) 

P 2 (4 1 1) X 2 
(. 05 ) ( 10 ) 

P 2 ( 5 lD X 2 
1 ) ( 10 ) 

P 2 ( 6 |Dx 2 

P 2 (7 1 1) X 2 

P 2 (8 1 1) T 2 
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The construction of the full set of mass functions 
results in the constant matrix A shown in Figure 7, and the 
functional matrix [si - B] in Figure 8. 

The method of solution is to invert [si - B] by the method 
of cofactors using the relation 



p . . (s) = 

j 



Isl - B 



M . . (s) 



where |sl - b| is the determinant of the matrix [si - B] , and 
Mj^(s) is the determinant derived from |sl - b| by deleting 
the jth row and ith column. 

The advantage of the cofactor method is that it permits 
the elements of [si - B] ^ to be computed individually saving 
time in cases where only a portion of the inverse is needed. 

The results are as follows: 



~ \ 1 T { “7ti 

= Lie i, 



s + 7 



Pl2(s) = 



2 _ 1 1 
(s+7) (s+5) s+5 s+7 



T / ~5t “7ti 

= Lie - e }, 






2 ^ _2 2 _ 

(s + 7 ) (s+6) s+6 s+7 



-6t 



= Ll2e“°^ - 2e“'^} 



-7t,, 



Pl4(s) = 



(s+7 ) (s + 6) s+6 s+7 



^ r -6t -7t, 

= Lie - e ; , 
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0 

0 

0 
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0 

3 
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0 
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1 
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2 
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4 

4 



Figure 7. Constant Matrix A 
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-2 
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-1 
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0 
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-1 
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-2 
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0 


s+4 
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-4 
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-4 
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s + 4 
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0 
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Figure 8. Functional Matrix [si - B] 
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(¥+7 )Ts+6) (s + 4 ) s (s+7 ) (s+5) (s+4) s (s + 7) (s+6) s 

. 2 , 4 1 

(s + 7 ) (s + 5) s (s+7) (s+4) s '(s+7 ) s 

1 2 ^ _3_ ^3 5 



s+7 s+6 s+5 s+4 



= L{l“2e + 3e ^^ + 3e - 5e 



As required ^ . (t) =1, t ^ 0 



j = l 

The interval probabilities (t) , . . . , p^^g(t)} are now 

known and the quality function can be computed for any 
functions f(t), and g(x) or d(i).|| 
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It should be noted that the example was particularly 
simplified by the fact that the matrix [si - B] had no non- 
zero entries below the diagonal. The reason for this is 
that all three components were considered to be nonrenewable. 
Thus when the unit left any state it could never return to 
it. And by arranging the states in order so that the earlier 
states corresponded to the higher rows, the triangular form 
was assured. The resulting denominator of all the expres- 
sions p^j (s) was of the form 

8 

|sl - B| = 7T (s+X ) 
i=l ^ 



which is already in the factored form necessary for the partial 
fraction expansion. 

If however at least one of the components had been renev;- 
able, there would- have been at least one positive entry below 
the diagonal and the denominator of the expressions p^^ (s) 
would have been of the form, 



8 

|sl - b| = 77“ (s+X . ) - K(s) 

i=l 



where K(s) is a polynomial function of s of degree less than 
eight. Therefore |sl - b| is a polynomial of degree eight 
whose roots and factors are unknown in general. This computa- 
tional problem is formidable since the degree of polynomials 
to be solved will equal N = 2^ (v;here n is number of system 
components) , and there is no exact solution for polynomial 
equations of degree greater than four. There are however 
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computer routines which perform approximation procedures, 
which may prove to be of value in employing this model. 

Finally, a rough measure of effectiveness may be computed 
from q(t) without inverting the Laplace transforms when 
f (t) =1, t ^ 0. Let the quality index be defined as follows; 




N 

j=i 



Thus Q is the expected aggregate value over an infinite time 
period. Of course, any realistically defined unit will 
reach a state with zero effectiveness in finite time, so 
that Q has a finite value. However this measure does not dif- 
ferentiate between units whose quality functions are thinly 
distributed over the real time line and those v.’hose effec- 
tiveness is heavily concentrated in short intervals of reason- 
able engagement duration. 
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V. CONCLUSION 



The purpose of this thesis has been to outline a proposed 
modeling viewpoint. It has not been advanced as a finished 
product. It is hoped that the concepts introduced will be 
considered as a worthwhile point of departure for additional 
research, which could result in establishment of a valuable 
decision tool for choosing among alternative unit structures. 

The most critical unsolved problem is the computational 
difficulty involved in solving the model. In addition to the 
inversion of Laplace transforms, the current solution neces- 
sitates the evaluation of large order determinants, which can 
become unwieldy even with electronic computers. Thus the 
current span of application is limited to small organizations. 
Research to extend this span should concentrate on finding 
simplified computational procedures for exact solution, or 
alternatively, determining efficient and accurate approximation 
methods. The use of the infinite series definition of the 
matrix exponential function as an approximation method was 
introduced in Section IV. Unfortunately, this expression 
requires sequential multiplication of large order matrices, an 
operation which is as computationally explosive as evaluation of 
determinants . 

Current research in reliability theory may offer a valuable 
approximation method, particularly since all of the random 
processes and measures of effectiveness were shown to be 
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expressible as linear functions of variables related to 
reliability of coherent systems. Esary and Proschan [Ref. 4] 
show the usefulness of reliability bounds and approximations 
based on the concepts of minimal cut sets and minimal path 
sets of components. Adaptation of this method to the avail- 
able type of data inputs may prove useful to the current model. 

A possible exact solution method using flow graph analysis 
is proposed by Howard [Ref. 6], but at present it appears to be 
subject to an even greater computational burden than the cur- 
rent method. 

In order to enhance the accuracy of the model it may be 
desirable to alter the assumption of exponentially distributed 
interevent times in favor of some other distribution. This 
would change the process to one of the semi-Markov variety. 

These processes are more difficult to solve in general, but 
offer greater flexibility in definition. Ross [Ref. 11] and 
Howard [Ref. 6] consider several solution methods. 

A highly desirable addition to the model would be a provi- 
sion for random threat. Rather than assuming constant Poisson 
rates of enemy fire, such a refinement would allov,-- rates to 
vary randomly, thus representing periods of non-engagement as 
well as change of situation during engagem.ents . The analogous 
reliability concepts of random wear and random shock have been 
studied by Gaver [Ref. 5] and Reynolds and Savage [Ref. 10]. 

Despite its capacity for refinement, the proposed model in 
present form does fill many of the needs for the study of struc- 
ture. It offers the decision-maker a model which uses his own 
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judgment, in an intuitively clear process of parameterization, 
in order to compute a functional measure of effectiveness 
based on fundamental functions which he, himself, defines. 
Indeed, with the exception of clearly stated assumptions made 
in building the m.odel, and the purely mathematical solution 
process, the model in use is completely established by the 
user, with whom ultimate decision responsibility rests. 
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